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ABSTRACT 

We have used archival far-uhraviolet spectra from observations made by the Space Telescope Imaging 
Spectrograph (STIS) of the Hubble Space Telescope and the Far Ultraviolet Spectroscopic Explorer 
(FUSE) to determine the column densities and rotational excitation temperatures for carbon monoxide 
and molecular hydrogen, respectively, along the lines of sight to 23 Galactic O and B stars. The 
sightlines have reddening values in the range E{B — V) = 0.07-0.62, thus sampling the diffuse to 
translucent interstellar medium. We find that the H2 column densities range from 5 X 10^^ - 8 x lO^o 
cm~^ and the CO from upper limits around 2 x 10^^ cm~^ to detections as high as 1.4 x 10^^ cm^'^. CO 
increases with increasing H2, roughly following a power law of factor ^ 2. The CO/H2 column density 
ratio is thus not constant, and ranges from 10"'' — 10~^, with a mean value of 3 x 10~^. The sample 
segregates into "diffuse" and "translucent" regimes, the former having a molecular fraction less than 
^ 0.25 and Ay /d < 1 mag kpc^^. The mean CO/H2 for these two regimes are 3.6 x 10"'' and 9.3x 10~^, 
respectively. These values are significantly lower than the canonical dark cloud value of 10~*. In six 
of the sighthnes, the isotopic variant ^^CO is observed, and the isotopic ratio we observe 50 — 70) is 
consistent with, if perhaps a little below, the average ^^C/^'^C for the interstellar medium at large. The 
average H2 rotational excitation temperature is 74 ± 24 K, in good agreement with previous studies, 
and the average CO temperature is 4.1 K, with some sightlines showing temperatures as high as 6.4 
K. The higher excitation CO is observed with higher column densities, consistent with the effects 
of photon trapping in clouds with densities in the 20-100 cm~'^ range. We discuss the implications 
for the structure of the diffuse/translucent regimes of the interstellar medium and the estimation of 
molecular mass in galaxies. 

Subject headings: ISM: abundances, ISM: clouds, ISM: lines and bands, ISM: molecules, ISM: structure 



1. INTRODUCTION 

Molecular hydrogen (H2) is the most abundant 
molecule in the interstellar medium (ISM), residing pri- 
marily in the large complexes of the dense molecular 
clouds that account for 10-20% of the mass in the inner 
disk of the Galaxy (jShull fc BeckwithlfTosl l. However, 
it is difficult to observe directly. H2 is a homonuclear 
molecule, with quadrupolar ground-state transitions that 
emit radiation only very weakly. Carbon monoxide (CO) , 
in contrast, has strong ground-state transitions that pro- 
duce readily observable emissions at radio wavelengths. 
In clouds with densities above the critical density for CO, 
H2 collisions dominate the ground-state excitation, and 
thus CO radio emission acts as a tracer of H2 in these 
regions. 

This relationship is often characterized by the conver- 
sion factor X — iV(H2)//cOi where Iqo is the integrated 
brightness temperature of the J = 1 — radio emission 
line at 2.6 mm. The value for the conversion factor is 
generally determined by one of the following techniques 
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(cf. lYoung fc ScovilldHMI) : correlation of the CO emis- 
sion with Ay in clouds determined by star counts, which 
is then correlated with H2 by the extrapolation of the 
Nh/Av from the diffuse ISM (Savage et al. 1977); an ex- 
citation analysis of ^^CO, assuming it is optically thin, 
i^CO is optically thick, and the ^^CO/^^CO is known; a 
virial analysis using the cloud sizes a nd linewidths; and 
comp arison with 7-ray emission (e.g. IStrong fc Mattoxl 
|1996[) . These techniques produce a value for X of about 
2 X 10^° cm~^ (K km s~^)~^, but examples of up to a 
factor of 10 higher can be found in the literature. 

It is generally recognized that CO emission can qualita- 
tively trace the distribution of H2 in the ISM, but because 
of difficulties in translating a measurement of /co(l ~ 0) 
into a column density of CO, X is not a good quantita- 
tive measure of 7V(CO)/iV(H2) (hereafter CO/H2). To 
relate X to a relative column density ratio requires the 
assumption of LTE and is typically done with the un- 
saturated ^^CO radio lines, so one must also assume ex- 
citation temperatures and abundance ratios for the two 
CO isotopes. This exercise gives a value of CO/H2 of 
about 10"** for the de nse molecular clouds, but with a 
large uncertainty (e.g. lDickman|[l978l ). 

The determination of CO/H2 is simplified by the di- 
rect observation of the two species in absorption. This 
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Fig. 1. — Sample CO absorption profiles for sightlines with increasing column density. The best fitted values for column density, rotational 
excitation temperature and Doppler broadening parameter are listed. 



has been done in the infrared, using vibrational tran- 
sitions, and values of 2 — 3 x 10~^ vifere found 
for the molecular clouds in NGC 2024 and NGC 2264 
(|LacY et al.l 11994) . With ultraviolet observations, this 
direct method can be extended to the diffuse and 
translucent phases of the ISM, where the molecules 
are more susceptible to the effects of photodissociation, 
and thus CO/H2 is expected to vary relative to the 
dense clouds. Models suggest that the photodissocia- 
tion of CO is sensitive, even more so than H2, to the 
strength of the interstellar ultraviolet radiation field, 
cloud geometry and the ultrav iolet abso rpti on and scat- 
tering properties of dust ([ van Dishoeck & Black' [19881 : 
iKopp. Roueff. fc Pineau des Forcts 2000 ) . Therefore the 
inter-relationship between CO and H2 may also be a good 
measure of the physical conditions and structure of the 
diffuse and translucent regimes of the ISM. Furthermore, 
there is evidence that a significant contribution to the 
large-scale Ga lactic CO emissio n is made by lower opti- 
cal depth gas (jPolk et al.lll988f l. 

H2 is best observed in the far-ultraviolet where there 
are extensive dipole-allowed electronic absorption band 
systems. The Far Ultraviolet Spectroscopic Explorer 
{FUSE) has proven to be an excellent instrument 
for studying H2 in absorption fRachford et al.l 120021 : 
iTumlinson et all 120021 : iGillmon et al.ll2006 l. aiidothers). 
CO also has an absorption band system throughout the 
ultraviolet, and a study of CO in the diffuse ISM us - 
ing Copernicus was perforr ned bv | F^derman et al.l (|1980D 
and later reanalyzed by iCrennv &: FedermanI (|2004D . 
These observations were limited to the C-X (1088 A) 



and E-X (1076 A) absorption bands toward 48 nearby 
bright stars with only about one third having CO, some 
of those being saturated. Further, the spectral resolu- 
tion {R ~ 20000) was not high enoug h to resolve the 
rotati onal substructure of the bands. iFederman et al.1 
(|1980[ ) concluded that measurements with a more sen- 
sitive, higher spectral-resolution instrument would allow 
for a better understanding of the relationship between 
H2 and CO across a wider range of environments. 

This can be achieved using the Space Telescope Imag- 
ing Spectrograph (STIS) on board the Hubble Space Tele- 
scope, which has access to the Fourth Positive (A^II - 
X^S+ v" — 0) band system. More than a dozen absorp- 
tion bands from the ground vibrational state of CO have 
been detected and, because of the large variance in os- 
cillator strengths, they probe a variety of optical depths, 
allowing for an accurate determination of column den- 
sity. Additionally, with the highest resolution grating 
modes (E140H) these bands can be resolved clearly into 
their constituent rotational structure. Because the rota- 
tional levels in the ground state of CO are closely spaced, 
the rotational excitation temperature of the molecule is 
e asily determined. 

IPan et al.l ((20051) used STIS Echelle and FUSE data 
to explore the CO and H2 absorptions in the Cepheus 
0B2 and 0B3 clusters. They find evidence for systematic 
variations in CO/II2 in these two different star-forming 
regions, which they say may indicate differences in star- 
formation histories. In this paper, we broaden the study 
of CO/H2 to the diffuse molecular regime of the ISM, 
rather than study isolated regions. We present an analy- 
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sis of 23 stars, which have been observed by both FUSE 
and the STIS E140II mode. The reddenings of these 
stars range from E{B — V)= 0.07 — 0.62, complementing 
the Pan et al. study well, whose sightlines range from 
E{B - V)^ 0.35 - 0.86. The Ha and CO absorptions 
have been measured and column densities and rotational 
temperatures determined. These data provide a survey 
of the CO/H2 relationship in the diffuse to translucent 
molecular regime of the ISM. 

2. DATA AND ANALYSIS 

The data presented here were retrieved from the Mul- 
timission Archive at STScI (MAST). The STIS observa- 
tions employed the E140H grating, providing the highest 
possible spectral resolution. The wavelength coverage of 
STIS allows for the observation of a number of absorp- 
tion bands of the CO A-X {v" — 0) band system. For 
some of the stars, though, the specific tilt of the Echelle 
grating allows for observation of only higher vibrational 
bands {v" > 7), because of the more blueward wave- 
length coverage. The FUSE data were processed with 
the CalFUSE pipeline, version 2.2. The individual chan- 
nels (e.g. LiF la, SiC la, etc.) were joined using an IDL 
shifting routine written to combine time-tagged or his- 
togrammed FUSE data. Additionally, we obtained STIS 
data corresponding to five sightlines from a, FUSE survey 
of translucent clouds (jRachford et al.ir2002l ). from whom 
we obtained the H2 column densities. 

2.1. Carbon Monoxide 

CO has an extensive electronic absorption band sys- 
tem, the Fourth Positive {A^H — X^'S'^) system, ranging 
from 1510 A to shorter wavelengths. Because the ener- 
gies of the rotational levels in the ground state of CO are 
closely spaced, the relative strengths of the individual 
ro- vibrational transitions are very sensitive to excitation 
temperature. Previous surveys, such as Federman et al. 
(1980) used Copernicus data, which were not able to re- 
solve the rotational structure but the high spectral reso- 
lution of the STIS E140H grating (i? = A/AA = 110, 000) 
allows these transitions to be well-resolved from each 
other and provides for an accurate determination of the 
rotational temperature in addition to column density. 

Model spectra of the ro-vibrational absorption were 
ge nerated using the wavele ngths and oscillator strengths 
of (jMorton fc Noreaul[l994D . The best fit model was cho- 
sen using a statistic from a grid corresponding to 
varying column density (N), rotational excitation tem- 
perature (Trot), and Doppler line broadening parameter, 
b. For most of the low-column sightlines, the bands (1- 
0) through (5-0) were fit simultaneously. For sightlines 
with only the shorter wavelength data, the bands start- 
ing with (7-0) and higher were fit. Figure [T] shows a sam- 
ple of spectra with five bands of the A-X system range in 
column densities with corresponding best fit models over- 
laid. When the lower vibrational bands were available, 
the detection limit was typically N ~ 2.5 x 10^^ cm~^. 
For the (7-0) band, it was ^ 4.0 x lO^^ cm^^ 
but four of the sightlines, CO absorption was detected, 
and in six of the higher column sightlines the isotopic 
variant ^^CO was also observed. 

In addition to the profile fitting, a curve-of-growth 
(COG) analysis was performed for each sightline. The 




1048 1050 1052 1054 1056 1058 1060 1062 
Wavelength (A) 



Fig. 2.— Sample H2 absorption (HD 93205) of the Lyman (4- 
0) band system, with best fitted absorption profile overplotted. 
Absorptions out of the J = — 4 are included in the fit; however, 
the J = and f lines contain the majority of the column density, 
and only these are considered in determining the total H2 column 
density. Other lines in the spectrum are from Ar I (1048. 2A), Fe II 
(1055. SA), and high J states of the Lyman (5-0) band. 

individual ro-vibrational absorption profiles of each avail- 
able band were fitted with Gaussians, their equivalent 
widths measured, and a COG constructed. A similar 
grid-search process as used in the profile fitting was em- 
ployed and the N, T^ot, and b determined by comparison 
to the theoretical COG for a single velocity component. 
For sightlines with column densities below lO^'' cm~^, 
the best fitted values agree with those determined from 
the profile fitting process within the error. 

For the higher column sightlines, unresolved velocity 
components may infiuence the results. The equivalent 
width of an absorption line can be higher for a given to- 
tal column density if the velocity structure contains un- 
resolved non-overlapping absorption lines. This effect is 
more pronounced as the lines become saturated, and thus 
the equivalent widths of the higher oscillator strength 
bands are inflated, producing a COG that mimics that 
of a single component with higher b. However, if enough 
unsaturated lines can be observed, then the column den- 
sity determination is still robust. The combination of 
the STIS signal-to-noise ratio and high spectral resolu- 
tion allow, in virtually all cases, for the observation of 
absorption lines that lie on the linear part of the COG. 
We noted that, for values of b above ^ 0.6 km s~^, sat- 
uration effects are mild to negligible for absorption lines 
with an NfX < lO^^ cm'^ A. 

Our initial fits of the higher column sightlines produced 
6-values in the range 1.0 — 3.0. These are higher than 
expected for molecular material (jPan et al.l (f2004) find 
b = 0.6 — 1.0 for CN, for example). We found that 
by limiting the fitting process to only the weaker ab- 
sorptions, we were able to reproduce the lower b values. 
In these cases, the fit was limited to those lines with 
W\/X ~ 5 — 9 X 10~®, depending on the data quality, to 
avoid saturation effects. The application of the equiva- 
lent width cut lowered the determined b value, and in- 
creased the determined N by typically 0.2 dex. The need 
for this approach was demonstrated best by our attempts 
to fit the HD 24534 (x Per) sighthne, the N{CO) of 
which was increased by 0.5 dex after the application of 
the equivalent width cut. Twelve bands of the A-X sys- 
tem were clearly observed, but we could get no consistent 
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Fig. 3. — CO (left) and H2 (right) column density versus reddening, E{B — V). 



fit from profile fitting. Also, the COG shows deviations 
from what would be expected for a single set of fit pa- 
rameters. This is most likely because of the presence of 
unresolved velocity components. 

The CO absorption profiles for x Per ha ve been fit be- 
fore u sing GHRS data of the A-X system (jKaczmarczykI 
120001) as well as STIS data of the i nter system bands 
( Sheffer. Federman. fc Lambert) I2002D . iKaczmarczvkl 
(2000) found that a two-component model, with one 
component having 85% of the column, was necessary to 
get a good profile fit of the data. He derives a column 
of N{CO) = (1 ± 0.2) X 10^^ cm-^ consist ent with 
that found by ISheffer. Federman. fc Lambert! (1.41 x 
10^^ cm-2) as well as this study (1.35 x 10^^ cm^^). 
In principle, the unsaturated inte rsystem bands, as 
ISheffer. Federman. fc Lambert! ([2002) used, could also be 
used for the highest column density sightlines, but there 
were no STIS E140H data available in the archive with 
the appropriate wavelength coverage for any of the other 
stars. 

2.2. Molecular Hydrogen 

H2 has an extensive band system in the FUSE wave- 
length range (912 < A < 1110 A at temperatures typ- 
ical of the diffuse ISM) arising from electronic transi- 
tions out of the ground state (X^S^) to the Lyman 

(B^T,^) and Werner (C^H^) levels. These absorption 
lines were first used by Copernicus to characterize the 
average properties of the molecular phase of the ISM 
(jSpitzer et all 11974!) . More recently, FUSE has been 
used to c ompile large samples of H2 absorption dat a for 
galactic ()Rach ford et a ll 120021: rGillmon et ani2006D and 
extragalactic flTumlinso n et al. sightlines. 

We have analyzed the H2 absorption lines in the sight- 
lines presented here with a least-squares fitting routine 
that compares the data to a model spectrum. The model 
spectrum was constructed with a power-law fit to the 
stellar continuum over the wavelength region of interest 
(w 1035 - 1140 A) with absorption compone nts derived 
using the H200IS optical depth templates of .McCandliss! 
(|2003! ). The continuum was placed interactively as in- 



dividual adjustments were necessary in almost all cases. 
The optical depth templates are calculated for each line 
using a Voigt profile, ensuring the correct line shape for 
different column density regimes. Figure [2] shows a sam- 
ple absorption spectrum of H2 and the best fitted model. 

An iterative approach is used for absorption lines from 
rotational states from J = to J = 4. A coarse grid of 
column densities (15.6 < log A^(H2( J)) < 23.0) is scanned 
for 6- values ranging from 2 to 11 km s~^. This process 
establishes the column to within 0.25 dex and 6-value 
to within ± 1 km s~^. A fine grid of column densi- 
ties (AA^ = 0.01 dex) is then scanned and the best fit 
is quoted to the nearest 0.05 dex. We feel that this is 
a conservative estimate of the accuracy of our H2 col- 
umn density determination. As a consistency check, we 
compare our re sults for two sightlines also studied by 
iRachford et al.l (^2002). The column densities measured 
for HD 185418 are identical to theirs within the errors. 
On the other hand, HD 102065 is discrepant in log A^(0) 
by 0.2 dex. 

This iterative approach proved repeatable for deter- 
mining column densities for the J = and 1 lines. The 
higher-lying rotational states (J = 2 - 4) were fit with 
the same procedure, initially scanning a grid of column 
densities from 10^'^ — 10^" cm~^. The quantitative results 
for these lines were less certain. Most of the lines arising 
from J = 3 and 4 fell on the "fiat" part of the curve of 
growth (as well as those from J = 2 in a few cases), and 
were subject to the degeneracy between column density 
and 6-value. Additionally, some of the lines from J — 2 
and 3 were blended, and we found that continuum place- 
ment was even more significant when dealing with these 
lines. To avoid these complications, we only present the 
columns derived for the J = and 1 levels as these states 
contain the majority of the molecular mass at diffuse ISM 
temperatures. 

3. RESULTS 

3.1. Column Densities and CO/H2 

Table 1 summarizes the derived column densities and 
rotational excitation temperatures for the CO and H2 
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along the line of sight to each star, as well as the spec- 
tral type and reddening. For all but three of the sight- 
lines in this study, we have obtained from the liter- 
ature values for the neutral hydrogen column density; 
the H I column densities for HD 27778, HD 102065 and 
HD 203532 were determined from Lya fits to the STIS 
data. The spectral type for HD 102065 given in the 
literature of B9IV is inconsistent with the presence of 
weak but notable C IV AA1548-1550 features typical of 
spectral types B1.5V or slightly later for high luminosity 
class (, Walborn et al.l [T995.1 . Consequently, wc adopt a 
B2V spectral type for this object ((B - V)o = -0-24), 
yielding an E{B - V) = 0.31. Table 2 hsts the col- 
umn densities of the first two rotational states of H2, 
the column density of H I and the molecular fraction, 
/ = 27V(H2)/(2iV(H2)+A^(H I)). 

With both molecules there is a trend toward higher col- 
umn with increased reddening (see Figure[3]), though nei- 
ther exhibits a tight correlation; for any given E{B — V), 
the column density of both molecules can vary by at 
least an order of magnitude. The scatter is most likely 
a product of varying environments. The color excess is 
a measure of the total amount of dust along the line 
of sight, but may not necessarily reflect the geometri- 
cal distribution of the gas and dust; a single translu- 
cent cloud could produce the same reddening as a se- 
ries of diffuse clouds - see, for example, the models of 
iKopp. Roueff. fc Pineau des ForetsI (|2000D . 

Figure [4] shows the correlation of CO with H2 . This 
is the tightest correlation seen in this study, roughly a 
power law relationship, i.e. N{CO) oc Af(H2)", with 
a w 2. At lower H2 columns, the CO can be difhcult 
to observe, whereas along the higher column sightlines, 
CO increases very quickly. This is most likely an indi- 
cation of the onset of CO self-shielding, but may also 
be attributed to increased shielding of CO from pho- 
todissociation by either H2 or dust. Overplotted are th e 
results from the survey by ICrennv fc FedermanI (|2004D , 
which used Cope rnicus data and the Cep OB data from 
iPan et all (j2005l ). who used STIS and FUSE data simi- 
larly to this work. Pan et al. note a systematic offset of 
the CO/H2 relationship between the Cep 0B2 and 0B3 
clusters, further indicating that differences in the physi- 
cal conditions of the clouds are playing a strong role in 
the scatter in this relationship. 

Contours show lines of constant CO/H2. The slope 
of the correlation is steeper than 1 and thus there is 
an increasing trend of the CO/H2 with increasing col- 
umn, with a steepening of the slope near a CO column 
of w 10^^ cm~^, where the self-shielding of CO becomes 
important. The value of CO/H2 reaches about 10~^ for 
the highest column sightlines, and if the general trend 
is extrapolated to high H2 column agrees well with the 
va lue for dense clou ds of about 3 x 10~^ as determined 
bv lLacv et aP (fl99l for NGC 2024 and NGC 2264. This 
trend is well supported theor e tically , such as by the mod- 
els of Ivan Dishocck fc Black (1988), which suggest that 
the column density ratio should increase from around 
10^'' in diffuse clouds to lO^'' in dense clouds. 

The left panel of Figure[5]shows CO /H2 as a function of 
reddening. Again the results of earlier studies are shown 
overplotted. There is a general trend toward higher ratio 
with increased reddening. An interesting standout from 
this trend are the Cep 0B3 data, which appear to have 
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Fig. 4. — Correlation plot of CO with H2 for the data presented 
in this study as well as previous studies. Contours of constant 
CO/H2 are overplotted. 

far smaller CO/H2 for the E{B — V) than expected. 

The middle panel of Figure [5] shows the ratio with 
respect to the molecular fraction. Below a molecular 
fraction of about 0.25 (the average value found in the 
ISavage et aLlll977l studv). the CO/H2 is about 3 x lO"'^ 
on average; we will consider this the "diffuse" regime. 
These lines of sight have a significant spread in redden- 
ing, covering 0.07 < E{B — V) < 0.46, and are not 
well differentiated from the high CO/H2 sightlines in the 
plot versus reddening. The fact that they separate well 
from the high ratio sightlines in this plot suggests that 
the molecular fraction represents a better measure of a 
cloud's physical condition than E{B — V). This is most 
likely because the amount of dust tends to follow the total 
hydrogen along the line of sight rather than the molecular 
component and the destruction of these molecules is more 
sensitive to line-shielding than dust-shielding. Above a 
molecular fraction of 0.25 there is marked increase in the 
ratio with increasing molecular fraction, with an average 
value of about 7 x 10"^. 

We see a similar behavior of the CO /H2 with the tra- 
ditional density measure Ay/d, where Ay is the magni- 
tudes of extinction in the V ba nd and d is the distance to 
the star. IValencic et al.l (|2004[ ) have compiled a large cat- 
alog of extinction properties to reddened Galactic O and 
B stars, and there is an overlap of 15 stars with this study. 
The right pan el of Figure[5lshows CO/H2 vers us this den- 
sity measure. iJenniskens fc Greenberd ()1993[ ) found that 
Ay/d — 0.9 mag kpc~^ appears to separate the diffuse 
from the dense lines of sight, and our data show an in- 
crease in the CO/H2 above Ay/d ~ 1; however, we would 
state that the division is between diffuse and translucent, 
as the sightlines in our study with Ay/d > 1 mag kpc~^ 
still have CO/H2 about an order of magnitude less than 
that for dense clouds. 

Further comparison to the IValencic et al.l catalog pro- 
duces no significant correlations of extinction curve pa- 
rameters to the measured values of N or T j-pt for either 
CO or H2, or for CO/H2. The study of iBurgh et all 
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Fig. 5.— 
(Av/d). 

(|2000f ) indicated a correlation between tlie CO column 
density normalized by E{B — V) and the strength of the 
far-UV rise; however, that study, using data from the 
International Ultraviolet Explorer satellite, covered col- 
umn densities ranging from A^(C0)=1G^'' — 10^^ cm~^, 
the correlation only becoming apparent at the highest 
columns and steepest UV curves. It is likely that in the 
environments of the clouds probed in this study, shielding 
by H2 and self-shielding are more important in regulat- 
ing the abundance of CO than dust-shielding, which may 
not play a significant role until A^(CO)> 10^^. 

3.2. Isotonic Fractionation 

For six of the sightlines, absorption lines from the iso- 
topic variant ^'^CO are observed. Table 3 lists the ^^CO 
and ^^CO column densities and rotational temperatures 
for these sightlines. The first column of ^'^CO column 
densities are those determined if it is assumed that the 
two isotopes have the same rotational temperature. Also 
listed are the results if the temperature is allowed to 
be a free parameter. The results give column densities 
within 0.2 dex of each other. The next section will dis- 
cuss the observed difference in rotational temperatures, 
when that parameter is left free to be fit. 

The isotopic ratio ^^CO/^^CO ranges from about SO- 
TO, with an average value of 57 ± 7 and we see no strong 
correlations with any other measurable quantity. The 
largest value is for y Per (HP 24534), whose ratio was 
determined by ISheffer. Federman. fc LambertI (|2002D to 
be 73 ± 12, in good agreement with the value found 
here; however, we do not see the enhanced fraction- 
ation measured alo ng such lines of sigh t as those to 
o Oph A. Y Qph (iFedermanet al.l l2003f). and C Oph 
(jLambert et al.lll994D . with values of 125 ± 23, 117 ± 35, 
a nd ^ 170, resp ectively. 

iWilson fc Ro od (1994) review the literature for deter- 
minations of the interstellar ^^C/^^C and adopt a value 
of 77 ± 7; however there is significant scatter amongst 
the values reported in the literature that seems to de- 
pend on the methods used, wh ich can include from mil - 
limeter emissions of C^^O (e.g. iLans-er fc Pe"^^^ (IT993h 
get ^"^G/^^G— 57 — 74) and nea r-infrared sp ectroscopy 
of CO vibrational bands (e.g. iGoto et alT ((2003.) get 
i2c/i3C= 86 - 137). One method that is not sensi- 



tive to processes of selective fractionation (described be- 
low) is determining t he ^^CH+Z^'^CH"'" rati o and the re- 
cent measurements of lCasassus et al.l()2005D agree with a 
value of ~78 but they interpret the scatter in their data 
(Icr = ±12.7) as a true measure of chemical heterogeneity 
in the local ISM. 

There are two main processes that can cause 
^^CO/^^CO to deviate from the average ISM value of 
the ^^C /^^C isotope ratio. T he first is isotopic charge ex- 
change (Watso n et al.lfl976l ). which occurs in gas where 
C+ is in abundance, as might be expected in translu- 
cent clouds, and enhances the ^'^CO because of its lower 
zero-point energy. The other proces s is selective isotopic 
photodissociation (jBallv fc Langed [1982) , which favors 
the more abundant, and thus more likely to self-shield, 

12C0. 

If isotopic exchange is more important than 
photodissociation we would expect ^^CO/^^CO = 
exp(-Ai;/fcTkin)x(i2c/i3C), where A£' is the zero 
point energy difference between the two isotopes 
(AE/k — 35 K). The average kinetic temperature for 
these six sightlines is 58 K, and therefore we could 
expect to see as low as half the average ISM isotopic 
ratio. However, it is unlikely that photodissociation does 
not play a role, and the models of Ivan Dishoeck fc Blackl 
(1988) suggest at best only a mild relative increase in 
the ^"^CO abundance for temperatures in the range seen 
in this sample. 

3.3. Rotational Excitation 

Both the FUSE and the STIS data are of high enough 
resolution to resolve the individual rotational transitions 
of H2 and CO respectively. This allows for the deter- 
mination of the rotational excitation temperature. For 
both molecules the relative column densities in the J = 
and J — 1 states are representative of the rotational ex- 
citation temperature of the gas, Tqi, following: 

iV(l)/A^(0) = gi/goeM-Eoi/kToi) 

where N{J) is the column density in the Jth rotational 
state, gj is the statistical weight, Eqi is the energy dif- 
ference between the J = and J = 1 states, and k is 
Boltzniann's constant. 
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Fig. 6. — CO rotational excitation temperature versus H2 ro- 
tational excitation temperature. The vertical dotted line is at the 
average value for the H2 temperature. Only sightlines with below 
average H2 temperature exhibit enhanced CO temperature. 

Transitions between the ground-state rotational states 
of H2 are not dipole-allowed, and have lifetimes long 
enough that at the densities typical of diffuse and translu- 
cent clouds collisions dominate over other processes in 
determining the relative populations of the low-lying J 
levels. Thus, Tqi for H2 traces the kinetic temperature 
of the gas. In CO, on the other hand, the rotational 
transitions are dipole-allowed and have short lifetimes 
such that the rotational excitation temperature will re- 
flect mostly the density in the cloud, not the kinetic tem- 
perature. 

The average H2 temperature in our sample is 74 ±24 K. 
This agrees very well with the value derived from Coper- 
nicus observations of 61 stars with H2 columns greater 
than 10^^ cm-2 of < Tqi >= 77 ± 17 K (ISavage et al 
119771) as well as the < Tqi >= 68±15 K of lRachford et"aLl 
()2002| ). The higher J levels can be populated through 
collisions, UV- and formation-pumping, and radiative 
cascade and thus their relative populations may not re- 
flect the kinetic temperature. The absorptions from these 
levels will not be considered in this study and should not 
greatly affect the column density determinations because 
the J — and 1 absorption accounts for the vast major- 
ity of the II2 column. 

There is no correlation of the rotational temperatures 
of either H2 or CO with reddening along the line of sight, 
nor was any seen in either the Copernicus study or the 
FUSE study of Rachford et al. The rotational temper- 
atures are also independent of total H2 column density. 
However, there are several interesting relationships with 
r otational temperature. 

iRachford et al.l (|2002t ) point out that for lines of sight 
dominated by translucent clouds, the molecular fraction 
should be large while the kinetic temperature is small. 
Though we do not see a strong correlation (i.e. we see 
several sightlines with low temperature and molecular 
fraction) all of the sightlines with above average H2 tem- 
perature have a molecular fraction below 0.2, except for 



HD 185418, which has a molecular fraction of 0.49. 

For CO column densities less than 10^^ cm~^, the av- 
erage rotational excitation temperature is 3.6 ± 0.5 K. 
For the higher column density sightlines, the average is 
5.2 ± 1.0 K. Also, the average CO rotational temperature 
is 4.5 ± 1.1 K for sightlines with below average H2 tem- 
perature, and 3.5 ± 0.5 K for those above (see Figure [5]). 
The increase of excitation with column density may be 
a result of "photon trapping" in clouds of higher optical 
depth. As the optical depth of the 2.6 mm J = 1 — ra- 
dio line increases, it becomes more likely that the photon 
will be reabsorbed by another CO molecule. This hinders 
the radiative cooling and the coUisional excitation/de- 
excitation process with the surrounding H2 molecules will 
equilibrate at a higher rotational temperature. The CO 
Trot is thus sensitive to the H2 space density, the H2 
kinetic temperature, and the line center optical depth, 
which itself depends on the column density and 6-value. 

Using an on-line statistical equilibrium radiative 

transfer code, available at 

http : //www . strw . leidenuniv . nl/$\sim$inoldatal 
(Schoicr ct al. 2005), we determined that the increase in 
temperature for the higher column density sightlines is 
consistent with the "photon trapping" effect in clouds 
with H2 space densities in the 20-100 cm~'^ range, 
assuming our average H2 temperature as the kinetic 
temperature and a 6- value of about 0.6 km s~^. This 
code assumes a single, homogeneous cloud, which may 
not necessarily represent a sightline with multiple 
velocity components. 

Another method for increasing the rotational 
excitation of CO is that of radiative pump- 

ing from nearby dense clou ds, as studied by 

IWannier. Penprase. fc AnderssonI (jl997[ ). This ef- 
fect will be proportional to the solid angle subtended by 
the nearby cloud as viewed from the CO along the line of 
sight, when the relative velocity between the absorbing 
and emitting clouds is zero. They point out that for 
sightlines with space densities in the few hundreds per 
cm^, and with kinetic temperatures in the range of the 
H2 temperatures we observe, fractions of '--^0.2 of An 
steradians is enough to account for the increased CO 
t emperatures observed. 

IWannier. Penprase. &: AnderssonI (|1997| ) point out 
that this effect may be differentiated from collisional ex- 
citation by observing other isotopes of CO, because the 
lower brightness temperature of the nearby clouds will 
reduce the rate of radiative pumping. Table 3 lists the 
^•^CO temperatures derived in this study, with the iso- 
topic temperature ratio T12/T13 listed in the final col- 
umn. For two or possibly three of the sightlines, there is 
an increased relative temperature for the more abundant 
isotope. These are also the sightlines with the higher 
molecular fractions, and are perhaps sightlines that pass 
through the outer portions of a denser molecular cloud, 
where a larger fraction of its sky may be covered by op- 
tically thick, radiating gas. However, we note that the 
"photon trapping" could also explain the T12/T13, be- 
cause the less optically-thick ^■^CO would not have its 
rotational temperature enhanced by the effect. 

4. DISCUSSION 

The traditional view of the ISM differentiates be- 
tween diffuse and translucent clouds, where the transi- 
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tion occurs at about Ay = 1. However, recent s tudies 
(jRachford et all 120021 : ISonnentrucker et alJ I2003D have 
suggested that there is not a clear distinction, and ar- 
gue that perhaps none of the sighthnes they see have 
truly "translucent" clouds, which models expect to show 
molecular fractions reaching very nearly 1. This situ- 
ation could arise from having multiple "diffuse" clouds 
along the line of sight, adding up to a given Ay- This 
would result in E{B — V) not necessarily being a good 
measure of the physical properties along a given line of 
sight, and our study confirms this. 

It is thus more important to measure parameters that 
are better representative of the physical state of the 
clouds. The molecular fraction gives one such mea - 
sure, but for reasons laid out in 'Rachf ord efall (120021 ) 
may not allow for the unambiguous separation of diffuse 
and translucent sightlines. However, the CO molecule is 
more sensitive to photodissociative processes, and can be 
used to probe the tr ansition region between diffuse and 
translucent clouds. iKopp. Roueff. fc Pineau des Foretd 
(|2000.) model the sensitivity of the abundance of CO to 
the effects of geometry, dust shielding and fragmenta- 
tion of the ISM. Their results suggest that if a sightline 
were simply a collection of small, diffuse clouds we would 
not expect to see an increased CO/H2 with total column 
density. The increase in CO/H2 that we observe with in- 
creased molecular fraction suggests that these sightlines 
do begin to sample the translucent regime. 

We believe that the transition from low to high CO /H2 
is similar to the one that hydrogen undergoes when the 
H2 column density gets high enough for self-shielding to 
take effect. The molecular fraction observed in a cloud 
will be dependent on t he balance of the fo rmation and 
the photodestruction. I Savage et al.l ()1977l ) found that 
the molecular fraction transitions from "low" (< 0.01) 
to "high" (> 0.01) values at about iV^ = 5 x lO^" cm'^, 
which corresponds to E{B — V) ^ 0.08 . This is also 
seen for high-latitude Galactic sightlines l|Gillmon et al.l 
|2006[) . The location of this transition is dependent on 
the physical conditions in the ISM, including the H2 
formation rate on grains, the intensity of the interstel- 
lar radiation field, and the degree of fragmentation of 
molecular clouds along the line of sight. The behavior of 
the Galactic sightlines is in marked coiitrast to that of 
the Magellanic Clouds ijTumlinson et al.]|2002f ). in which 
this transition occurs at much larger columns, perhaps 
because of the lower H2 formation rates and higher disso- 
ciating radiation fields. All of the sightlines in this study 
have total hydrogen columns above this transition value, 
and should be considered part of the "high" molecular 
fraction regime. 

The values of CO/H2 that we measure for the dif- 
fuse and translucent regimes range from 10"'^ - 10"^ 
in contrast with the canonical dense cloud value of 
10~'^. This discrepan cy is consistent with the st udies of 
Ide Vries et al] (|198/D and iMagnani et al.l (|1998D . which 
find variations in the X-factor as determined from radio 
observations of high-latitude translucent clouds. Both 
studies interpret their results as variations of the CO 
abundance, which could have a deleterious effect on the 
process of determining the masses of individual clouds. 
It could also impact the "weighing" of galaxies if the dif- 
fuse or translucent ISM contribute si gnificantly to the 
CO radio emission. I Polk et al] (|1988l ) found, based on 



the ratios of integrated -'^^CO and ^^CO emission, that a 
significant contribution to the total ^^CO emission from 
a galaxy could be arising from material of moderate op- 
tical depth. In that case, applying a standard X-factor, 
derived from the dense clouds, would result in the under- 
estimation of the total mass because of the lower CO/H2 
in the more diffuse gas. 

This difference could be particularly important 
for clouds subjected to intense UV radiation fields. 
lYao et al.l (|2003) found in their study of CO emission 
in starburst galaxies an X-factor significantly lower than 
the standard Galactic value. Further, they found, by 
virial analysis, that the clouds in these galaxies were not 
gravitationally bound unless the CO/H2 was 9-90 times 
lower than expected and suggest that the CO emission 
arises from nonviri alized warm and diffuse gas clouds. 
On the other hand, iRosolowskv et al.l (|2003l ) find no ev- 
idence in M33 for a diffuse molecular component traced 
by CO emission, concluding that most of the CO flux 
resides in the giant molecular clouds. 

The direct comparison of the column densities deter- 
mined from UV absorption line spectroscopy of CO and 
H2 provides a technique for probing the structure and 
chemical balance of the ISM free from many of the bi- 
ases and assumptions inherent in the traditional X-factor 
method. If the CO J=l-0 radio emission along these lines 
of sight could be measured then the relationship between 
CO/H2 and X-factor could be checked directly. It would 
also be useful to use absorption spectroscopy to explore 
the sightlines with iV(H2)= 10^^ — 10^^ cm~^, the regime 
in which we expect to see the transition from translucent 
to dense clouds. The power law increase in CO versus H2 
seen in the diffuse/translucent regime should turn over 
as Ay rises above ~ 3 and the CO/H2 "saturates" at the 
canonical dense cloud value of ~ 10~^. 

Furthermore, extragalactic environments could be ex- 
plored, such as the Magellanic Clouds; however, we do 
expect that with the currently available archival data 
CO will be diff i cult t o detect there. The study of 
iTumlinson et al.l (|2002f l found that although H2 column 
densities as high as around 10^° cm~^ have been de- 
tected, the typical columns found in the Clouds lie in 
the 10^^ - 10^^ cm-2 range. If the CO/H2 ratio were to 
follow the trend seen in the Galactic sample, the CO col- 
umn densities would be below the STIS detection thresh- 
old. We would expect, though, that the CO/H2 would be 
even lower because of the lower mctallicity and higher UV 
radiation fi elds seen in the Clouds. But, because mod- 
els predict ()Bell et al.ll2"00l e.g.) the CO/H2 to change 
with environment, if the appropriate data were taken, 
the technique could be used to explore the varying con- 
ditions in the diffuse/translucent ISM in other galaxies 
as well as our own. 
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TABLE 1 

Molecular Column Densities and Rotational Temperatures 









CO 


Ha 


CO/H2 








log JV 


^ rot 


b 


log iV" T,,t 


Star 


Sp. Type 




(cm^^) 


(K) 


(km s^'-) 


(cm-2) (K) 


(xlO-6) 



HD 24534*' 


BOVe 


0.45 


16.13±0.20 


5.3±0.6 


0.4±0.2 


20.92 


57±6 


16.1 ±8.3 


-LIJ--' 1 1 ) 


B3V 


0.38 


1 fi n^+n 1 4 


6.1it0.5 


7+0.3 

' -0.1 


20.79 


56±6 


18.0 ±7.1 




07V 


0.27 


< 13.60 






19.99 


61i7 


< 0.40 


HD 93205 


03V 


0.37 


13^23^10.06 


3.4±0.6 


4.3±1.0 


19.84 


105±21 


~" 0.24±0.07 


HD 93222 


OTHI 


0.40 


13.36±0.20 


2.8±1.0 


0.7±0.6 


19.80 


77±11 


0.36±0.19 


HD 93840 


Bl.OIb 


0.14 


13.38±0.08 


3.6±0.6 


-, -,+0.8 

—0 5 


19.55 


47±4 


0.68±0.21 


HD 102065'= 


B2V 


0.31 


13.62±0.12 


4.0±1.0 


1.7±i.2 


20.63 


59±7 


0.10±0.04 


HD 103779 


B0.5IH 


0.21 


<12.35 


- 


- 


19.90 


86±14 


<0.03 


HD 104705 


BOIb 


0.26 


12?98±0.16 


3.4±0.8 


1 0+i'2 


19.98 


92±16 


0.10±0.04 


HD 116852 


09IV 


0.22 


13.28±0.04 


3.0±0.3 


4+0-5 


19.85 


70±9 


0.27±0.07 


HD 121968 


BIV 


0.07 


<12.30 






18.70 


38±3 


< 0.40 


HD 152723 


06.5in 


0.46 


13?88±0.15 


4.0±0.3 


1.0±0.2 


20.33 


73±10 


0.36±0.15 


HD 163758 


06.5Iaf 


0.35 


13.42±0.05 


4.5±0.5 


1.3±0.5 


19.85 


86±14 


0.37±0.10 


HD 177989 


BOIII 


0.25 


14.62±0.17 


3.7±0.4 


0.7±0.1 


20.23 


52±5 


2.4 ±1.1 


HD 185418 


B0.5V 


0.50 


14.82±0.20 


3.2±1.1 




20.79 


105±21 


1.1 ±0.6 


HD 201345 


09p 


0.32 


<12.40 






19.43 


147±41 


<0.09 


HD 203532 


B3IV 


0.32 


15?70±0.17 


5.5±0.3 


0.7±0.1 


20.71 


49±5 


9.8 ±4.5 


HD 206267'' 


06.5V 


0.52 


16.11±0.13 


6.4±0.5 


Q 7+0.8 
-0.1 


20.86 


64±8 


17.9 ±6.8 


HD 207198'' 


09.5Ib 


0.62 


15.53±0.20 


3.7±1.2 


1.0±0.5 


20.83 


66±8 


5.0 ±2.6 


HD 210839'' 


Oeinfp 


0.56 


15.41±0.08 


4.4±0.3 


9+°-^ 


20.84 


72±10 


3.7 ±1.1 


HD 218915 


09.5Iab 


0.29 


13.64±0.13 


3.9±0.2 


1.6±0.6 


20.15 


86±14 


0.31±0.12 


HD 303308 


Bl.OIII 


0.30 


13.65±0.06 


3.7±0.4 


3.1±0.5 


20.35 


86±14 


0.20±0.05 


CPD -59 2603 


07V 


0.46 


14.15±0.10 


3.2±0.5 




20.15 


77±11 


1.00±0.33 



Errors are ±0.10*" H2 values from Rachford et al. (2002)'= Sp. Type and E{B - V) determined this study 
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TABLE 2 

Molecular and Atomic Hydrogen Column Densities 





H2 Column Densities (cm '^)° 


H I 




Molecular 


Star 


log Ar(0) 


log Ar(l) 


log N (cm~2) 


Ref. 


Fraction 


HD 24534° 


20.76 


20.42 


20.73 


1 


0.76 


HD 27778'' 


20.64 


20.27 


21.20 


4 


0.44 


HD 91824 


19.80 


19.55 


21.15 


2 


0.12 


HD 93205 


19.40 


19.65 


21.33 


1 


0.06 


HD 93222 


19.50 


19.50 


21.54 


2 


0.04 


HD 93840 


19.45 


18.85 


21.04 


1 


0.06 


rlD lUzUoo 


20.45 


20.15 


21.90 


4 


0.10 




19.55 


19.65 


91 1 


1 
1 


1 fl 




19.60 


19.75 


21.11 




0.13 




19.60 


19.50 






14 


HD 1 91 QfiS 


10. DO 


i / . I U 


20.71 




0.02 


HD 152723 


20.05 


20.00 


21.43 


1 


0.14 


HD 163758 


19.50 


19.60 


21.23 


1 


0.08 


HD 177989 


20.10 


19.65 


20.95 


1 


0.28 


HD 185418 


20.35 


20.60 


21.11 


2 


0.49 


HD 201345 


18.85 


19.30 


20.88 


1 


0.07 


HD 203532 


20.60 


20.05 


21.50 


4 


0.24 


HD 206267'' 


20.64 


20.45 


21.30 


3 


0.42 


HD 207198'' 


20.61 


20.44 


21.34 


1 


0.38 


HD 210839'' 


20.57 


20.50 


21.15 


1 


0.53 


HD 218915 


19.80 


19.90 


21.11 


1 


0.18 


HD 303308 


20.00 


20.10 


21.45 


1 


0.14 


CPD -59 2603 


19.85 


19.85 


21.46 


1 


0.09 



References. — (1) IDiplas fc Savaed l lT99i) . (2) IFitzpatrick fc Massal | |T990|) . (3) IRachford et al.l l l2002t) . (4) De- 
termined this study 
Uncertainties for N{0,1) ^ ±0.1*' H2 values from Rachford et al. (2002) 



TABLE 3 

CO IsoTOPic Column Densities and Rotational Temperatures 





1 i2co 


i^CO 








Log N 


Trot 


Log 


Log TV 


Trot 






Star 


(cm~^) 


(K) 


(cm~^) 


(cm-2) 


(K) 




712/T13 


HD 24534 


16.13±0.20 


5.3±0.6 


14.33 


14.30±0.12 


4.1±0.6 


68±31 


1.3±0.3 


HD 27778 


16.05±0.13 


6.1±0.5 


14.26 


14.28±0.08 


3.9±0.6 


59±14 


1.6±0.3 


HD 177989 


14.62±0.17 


3.7±0.2 


12.90 


12.82dz0.08 


3.9±0.5 


63±25 


1.0±0.2 


HD 203532 


15.70±0.17 


5.5±0.5 


13.95 


13.97±0.20 


4.8±0.6 


54±21 


1.1±0.2 


HD 206267 


16.11±0.17 


6.4±0.6 


14.44 


14.42±0.08 


6.4±0.8 


49±15 


1.0±0.2 


HD 210839 


15.41±0.04 


4.5±0.2 


13.75 


13.70±0.10 


3.5±0.8 


51±9 


1.3±0.3 



Column density assuming same Trot as for CO 



